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We study an exchange coupled system of itinerant electrons and localized fermion pairs resulting 
in a resonant pairing formation. This system inherently contains resonating fermion pairs on bonds 
which lead to a superconducting phase provided that long range phase coherence between their 
constituents can be established. The prerequisite is that the resonating fermion pairs can become 
itinerant. This is rendered possible through the emergence of two kinds of bond-fermions: individual 
and composite fermions made of one individual electron attached to a bound pair on a bond. If 
the strength of the exchange coupling exceeds a certain value, the superconducting ground state 
undergoes a quantum phase transition into an insulating pair-bond liquid state. The gap of the 
superfluid phase thereby goes over continuously into a charge gap of the insulator. The change-over 
from the superconducting to the insulating phase is accompanied by a corresponding qualitative 
modification of the dispersion of the two kinds of fermionic excitations. Using a bond operator 
formalism, we derive the phase diagram of such a scenario together with the elementary excitations 
characterizing the various phases as a function of the exchange coupling and the temperature. 

PACS numbers: 74.20.-z, 74.20.Mn, 03.75.-b 



I. INTRODUCTION 

The evolution of pairing correlations and, the related 
to it, onset of phase coherence in low dimensional sys- 
tems is at the center of intense theoretical investigations. 1 
This activity concerns systems such as: (i) high tem- 
perature cuprate superconductors with their spin/charge 
pseudogap phenomenon^ (ii) low temperature super- 
conducting materials which can be driven towards in- 
sulating or metallic phases via some extrinsic/intrinsic 
mechanisms j^^si (iii) ultra-cold gases of fermionic atoms 
in the cross-over regime between a BCS state and 
a Bose-Einstein condensate, controlled by a Feshbach 
resonance^. 

In this paper, we shall investigate systems where 
bosonic resonant pairs form in an ensemble of itinerant 
uncorrelated fermions. The dynamics of such a boson- 
fcrmion exchange coupled system is characterized by two 
competing processes: A) a local exchange between a lo- 
calized bound pair of fermions and a pair of itinerant 
uncorrelated fermions, B) a non-local single particle hop- 
ping of the itinerant fermions between nearest neighbor 
sites. 

The exchange between the localized bound pairs (act- 
ing as hardcore bosons) and the fermionic itinerant parti- 
cles creates a local quantum superposition with bonding 
and antibonding resonant pair configurations given by 
|2,i)± = -fcipf ± r+)|0). pf and r+ = c\^l denote re- 
spectively the creation operators of the two constituents 
of those local bound pair states, i.e. for hardcore bosons 
and pairs of itinerant fermions on site i. In the atomic 
limit, these states are separated by an energy difference 
equal to twice the boson-fermion exchange coupling, with 
the bonding state being energetically favorable. 

The local quantum structure of the bonding config- 



uration, responsible for inducing local resonant pairing 
among the fermions, has at the same time a hindering ef- 
fect in establishing a long range superconducting spatial 
correlations of such resonating pairs. Spatial correlations 
between the bonding/antibonding pairs are built up via 
other bosonic and fermionic bond configurations which 
form the local Hilbert space which has a very rich struc- 
ture (see Fig. la). There are two fermion- like configura- 
tions which result from attaching to the vacuum state one 
individual fermion (configuration "1"), respectively such 
an individual fermion together with a localized bosonic 
bound fermion pair (configuration "3"). Moreover, in 
this bond Hilbert space we have completely empty bonds 
- bond-holes "0" - and completely filled ones - double- 
pair-bond states "4". The fermionic objects are itiner- 
ant. Their mobility is due to the creation/annihilation of 
bonding (respectively anti-bonding) and bond-holes (re- 
spectively double-paired bonds). 

Our main aim in this study is to show that, due to 
the internal degrees of freedom of the two components 
(localized hardcore bosons and itinerant fermions), there 
are two distinct energy scales in this problem: one con- 
trolled by the dynamics of the bonding pairs and one by 
the dynamics of the bond-holes and double-paired bonds 
(see Fig. lb). Depending on the ratio between the ex- 
change coupling strength and the hopping amplitude, the 
dynamics of the two kinds of fermionic excitations can 
lead to either i) their pairing up in a superconducting 
state where simultaneously a coherent state with bonding 
pairs, bond-holes and double-paired bonds occurs, or ii) 
an insulating state with coherent bonding pairs and with 
zero amplitude of the bond-holes/double-paired bonds 
(see Fig. lb). 

The presence of bonding pairs does not guarantee by 
itself the occurrence of long range superconducting cor- 
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FIG. 1: The a) panel provides the description of the possible 
configurations on the single bond (hexagon). The bold arrow 
is associated with the pseudospin for the original localized 
bound pairs. The b) panel indicates the possible mechanisms 
of dissociation starting from a representative string of three 
neighbor bonding bosons. The successive two steps for that 
are: i) the formation of two kinds of pairs, 1 involving single 
individual fermions on a bond and the other made up with 
charges 3 due to individual fermions attached with a bound 
pair on a bond and subsequently ii) the bond fermions pairing 
up and resonating with the bond-holes and the double-pairs- 
bonds. 



relations between their paired constituents. They have 
to dissociate in order to induce a superconducting state. 
This process of dissociation involves two possible chan- 
nels for the dynamics of this resonant pairing scenario 
(see FigQJj): i) the single fermions and the compos- 
ite ones (constituted of a fermion-boson pair) hybridize 
and their itinerancy sustains a coherent liquid of bond- 
ing pairs, ii) the two kinds of fermions pair up and res- 
onate with the bond-holes and double-paired bonds and 
via that result in a superconducting state. 

Within this context, our target is the following: I) de- 
termine the evolution of the superconducting state into 
an insulating state, both quantum and thermally driven, 
II) analyze the nature of the excitation spectrum as a 
function of the exchange coupling. 

The outline of the paper is the following: in Sec II we 
present the scenario for resonant pairing, on the basis of 
three examples actively studied in the literature. In Sec. 
Ill we sketch the bond operator formalism and adapt it 



to the boson-fermion problem. Sec. IV is devoted to the 
derivation of the phase diagram as a function of exchange 
coupling and temperature features. In Sec. V, we present 
the evolution of the excitation spectrum as one tunes the 
ground state from a superconductor to an insulator. The 
last section VI is reserved for the summary, conclusions 
and outlook. 



II. THE BOSON-FERMION SCENARIO 

There is a variety of different physical systems where 
fairly localized bound states are quasi degenerate with 
itinerant states of their constituents. They can be para- 
phrased in terms of a two-component system composed 
of localized bosonic pair states, itinerant fermionic quasi- 
particles and a local exchange interaction between the 
two. Such a scenario is described by the so called boson- 
fermion model (BFM) Hamiltonian 



H = (D - M ) £ clc ia + (A B - 2 M ) 



- E *« (4<*r + H.c.) + g J2 {Ptn + P7r+) (1) 

where g is the strength of the exchange interaction and t 
the hopping integral for the itinerant fermions, which is 
assumed here to be different from zero only for nearest 
neighbor sites. The band half- width (which will serve 
as energy unit) is D = zt, z being the coordination 
number of the underlying lattice. The energy of the 
bound fcrmion pairs is denoted by Ab- The number 
of the ensemble of bosons and fermions being conserved, 
ntot = fiFf + up i + 2ns, implies a common chemical 
potential fi for both subsystems, rig and np a are the 
occupation numbers per site of the hard core-bosons and 
of the fermions with spin a ="f , J,. In the present study, 
we restrict ourselves to the analysis for the fully sym- 
metric half- filled band case ( Ab = 0), which in such a 
two-component system means ntot = 2. The annihilation 
(creation) operators for the fermions with a spin a at a 

certain site i are given by c^J , those for the hard-core 
bosonic fermionic bound pairs by the pseudospin- 1 pf 
and similarly those for the pairs of itinerant uncorrelated 
fermions by t+ = c| T c^, t~ = c a c lT . 

This model was introduced originally by one of us (JR) 
in an attempt to capture the salient features of polaronic 
systems in the intermediate coupling regime between adi- 
abaticity and anti-adiabaticity, but has turned out sub- 
sequently to be of much more general relevance and ap- 
plicability. Among others, such a scenario has led to the 
prediction of the charge pseudogap features in the high 
T c cuprate superconductors^, without having to invoke 
any particular microscopic mechanism for that. We shall 
discuss briefly three representative examples where such 
a scenario seems to be relevant. 

(i) In a system with strong local electron-lattice cou- 
pling we have the formation of small polaronic charge 
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carriers, which in general will exist in form of localized 
bipolarons. If their binding energy exceeds the band- 
width of the itinerant electrons in an undeformed lat- 
tice, this can give rise to bipolaronic superconductivity 

- an extremely fragile state of a phase fluctuation con- 
trolled condensate of bosonic tightly bound electron pairs 

- whose existence in real materials has yet to be experi- 
mentally verified. However, if the binding energy of such 
localized bipolarons is such that it overlaps with the con- 
tinuum of the bare electron band states, this will result 
in an exchange between such bound pairs and pairs of 
uncorrelated itinerant electrons close to the chemical po- 
tential. The microscopic mechanism for that exchange 
arises from large quantum fluctuations of the lattice dis- 
placements in the immediate vicinity of the charge car- 
riers. Such local lattice displacements fluctuate between 
essentially undeformed and much deformed lattice en- 
vironments which, as a result, periodically capture and 
release electron pairs on small local cluster (acting as ef- 
fective sites on a lattice) made up of atoms and their 
associated ligand environements^. 

(ii) In a system with strong local electronic correla- 
tions, as described by the Hubbard model close to the 
half-filled band case, hole pairing occurs on small pla- 
quette clusters^ This, coexisting with triplet spin pairs, 
destroys the antiferromagnetic long range order, giving 
rise to a spin liquid state made out of singlet hole pairs, 
which conceivably can condense into a superfluid state. 10 
One of the possible mechanism for that is the hopping of 
the singlet hole pairs between neighboring clusters. A 
more efficient way to get this condition is by exchanging 
the local bound hole pairs on a plaquette with single holes 
on neighboring plaquettcs. This mechanism results in an 
effective exchange coupling between bosonic bound hole 
pairs and pairs of uncorrelated fermionic holes on adja- 
cent plaquettes, which at the end assure the itinerancy 
in the system. Such a scenario has been analyzed on the 
basis of a kind of real space renormalization group tech- 
nique, where the initial Hubbard model can be rephrased 
in terms of an effective BFM, albeit with additional terms 
carrying the information of the underlying antiferromag- 
netic short range interaction— 

(iii) In a system with an optically trapped gas of 
ultracold fermionic atoms (studied in connection with 
the cross-over between a BCS state of weakly inter- 
acting fermionic atoms and a Bose-Einstein condensa- 
tion of tightly bound states of them) one can monitor 
the strength of the interaction between those fermions 

- sweeping it from a repulsive to an attractive interac- 
tion - via a so called Feshbach resonance mechanismiSii^ 

- under the effect of a magnetic field. This mechanism 
is based on binary collision processes in which the inter- 
atomic interaction depends on the electronic spin con- 
figuration of the pair. For an overall electronic spin 
triplet state (avoiding the Coulomb repulsion) this leads 
to a weakly bound state, while a singlet configuration 
leads to scattering states. Two incident atoms in a sin- 
glet configuration can enter into a resonance with such a 



weakly bound triplet configuration molecular state when 
their respective energies are close to each other. This 
is achieved by flipping the electron spin on one of the 
two atoms via a hyperfine interaction, thus acquiring the 
necessary triplet configuration to bind them momentar- 
ily into a pair. By applying a magnetic field during such 
binary scattering processes one can change the relative 
position of the energy levels of the two electron spin con- 
figuration and thus enter in a resonance regime where 
they are quasi degenerate. Such a situation can be then 
described by a phenomenological model such as the BFM. 

III. BOND OPERATOR REPRESENTATION 

To study the scenario for resonant pairing, sketched 
above, by explicitly taking into account the interplay be- 
tween the bosonic bonding pairs and the processes link- 
ing them to the single particle fermionic entities, we 
make use of the bond operator formalism and adapt 
it to the present BFM system. This approach turns 
out to be particularly appropriate in treating situa- 
tions where there is a natural pairing in form of dimers 
in the ground state, which is either imposed by the 
Hamiltonian (as in our case) or by a spontaneous sym- 
metry breaking. The bond operator theory has been 
successfully designed and applied in different contexts 
such as for antiferromagnetsi^i^, spin-ladder—, doped 
antiferromagnets: 1 '', bilayer quantum Hall"®, and Kondo 
lattice systems^. 

Introducing the bond operator formalism for the BFM, 
we start from a bond on each lattice site being made up 
of fermionic and bosonic configurations and then express 
the original fermionic and bosonic operators c^J , pf 1 , rf' 
in terms of the new basis. Concerning a single local bond, 
one has eight possible configurations to start with. Simi- 
lar to the case of a spin insulator, where one is introduc- 
ing the singlet and triplet boson operator, here we have 
four boson operators that describe two-particle configu- 
rations in form of bonding-bonds and antibonding-bonds 
and zero- as well as four-particle configurations in form 
of hole-bonds and double-paired-bonds. The two-particle 
bonding and antibonding configurations can be expressed 
in terms of pseudospin-i operators r^~, pf for the uncor- 
related fermions and the tightly bound ones. The zero- 
and four-particle configurations are described by |0) and 
Pi T i |0) respectively. We thus have to consider four bond 
boson operators 6^) ; a (t) j j (t) i aiK } ^(t) ( see pjg^ ^) with 

b f \v) = -±=[p + -r + }\0) 



P\v) = p+r+10) 

<ft\v) = |0>. (2) 

|0) = |0/) ® I Ob) denotes the boson-fermion vacuum, 
i.e., the state with no fermions (c a |0) = 0) and no 
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bosons present (/o~|0) =0). \v) designates a correspond- 
ing new vacuum state in which neither bond-bosons nor 
bond-fermions are present. The remaining configura- 
tions are four bond-fcrmion operators describing individ- 
ual fermionic states and states where such fermions are 
attached to a boson on a particular bond: 

h\\v) = 4|o) 

sUv) = P+410). (3) 

The operators , s*^) obey the canonical fermion com- 
mutation relations, while b"', a"' have boson 
commutation rules. In this representation, the total num- 
ber of states in the Hilbert space of these four bond 
bosons and four bond fermions is much larger than the 
eight configurations allowed by the BFM. To limit the 
kinematics to the physical subspace one has to impose 
the closure relation on the bond given by the constraint: 



IV. MEAN-FIELD PHASE DIAGRAM 

Let us now apply the bond operator method to study 
the spectral properties of the BFM, in conjunction with 
its various possible phases. We proceed by looking for 
solutions where the bosons associated with the bonding- 
bonds few, antibonding-bonds the bond holes 
and double-paired-bonds /(t) are condensed. This means 
that in the Hamiltonian, Eq^ after having transformed 
the original fermion and boson operators into the new 
bond operators, one is decoupling the quartic terms in 
a way such that the expectation value of the linear and 
bilinear bond boson operators have a non-zero amplitude, 
(o) = o, with o = 6, a, /, d. The coherent state which then 
emerges in the most general case has the following local 
structure: 

\ipc) = {u v + Ufc6 t + u CI a t + Uddfr + UfP)\v), (9) 



Pf + s Ua + tfb + ata + h\h a + dU = 1 . (4) 

a a 

In the subspace limited by the relation (Eq. 0J, one then 
derives the expressions for the original fermion and hard- 
core bosonic operators in terms of bond operators. 

Let us first consider the fermion creation operator with 
spin a. Making use of the closure relation (Eq. and 
considering all possible transitions induced by the single 
particle fermionic operator between the different bond 
quantum configurations, we obtain the following expres- 
sion: 

4 = Pa fU- a + -j=st {b + a) + ^= (at _ 6 t) h _ a + h l d , 

(5) 

where p a = (+1, —1) for a = (f , J.). In a similar way one 
obtains for the creation operator of a hard core boson: 

p+ = -L/t( a + £ s\K + ^ + at)d. (6) 

This, subsequently leads to the following expressions for 
the density operators for the bosons and fermions: 

p + p~ = ^a + b^b + a% + b^a) + pf + J2sUa 

(7 

4c = ^a + tfb-Jb-tfa)+pf+(sis„ + hih a ). 

(7) 

Finally, the expression for the pair exchange on the bond 
between the initial fermions and bosons reduces in the 
new representation to just the difference between the 
number of bonding and antibonding bosons, that is: 

p + r~ +T + p~ =a i a-bh. (8) 



where the finite amplitudes of the various coefficients 
u are related to the condensate of the corresponding 
bosonic degrees of freedom. Let us, at this point, make 
some general observations related to the structure of the 
Hamiltonian in such a bond operator representation. The 
6-bond bosons are always condensed at zero temperature 
(b 0), due to the local exchange and the hopping in- 
duced processes related to the creation of d-bond bosons 
and s, h bond fermions. Moreover, the d-bond bosons 
condense only if simultaneously the /-bond bosons con- 
dense and in our present analysis concentrating on the 
particle-hole symmetric half-filled band case, we have 
(d = / 7^ 0). In the present mean field approach the 
a-bond bosons have gapped excitations due to the ex- 
change splitting with respect to the 6-bond bosons and 
hence do not contribute to a coherent state in the lowest 
energy configuration. 

Let us now examine the possible solutions at zero tem- 
perature as a function of the exchange coupling g. Our 
main objective is to show how, above a critical value of 
the ratio g/D, the system passes from a superconducting 
state with simultaneous coherence of 6-bond bosons and 
the d-bond as well as /-bond bosons, to an insulating 
state where only the 6-bond bosons are condensed. 

Before going into the details about such a calculation 
and the results, let us point out certain aspects of the 
emerging bond boson- fermion dynamics which can be en- 
visaged on simple and quite general physical grounds. As 
mentioned before, the 6-bond bosons are always in a co- 
herent state at zero temperature, due the local exchange 
and the quantum processes of double creation and anni- 
hilation of nearest neighbor 6-bond bosons. Hence, the 
U(l) symmetry associated to the number conservation 
of the 6-bond bosons is generally broken in the entire 
regime of couplings g. In the weak coupling regime, b- 
bond bosons have a large dispersion which, together with 
coherent d-bond and /-bond bosons, leads to a super- 
conducting state. In the intermediate/strong coupling 
regime, the 6 bosons become more populated, which is 
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caused by an increased local binding energy. This in- 
crease in the b,d boson occupation leads to a decrease 
of the s- and /i-bond fermion densities, which ultimately 
causes a reduction of the b boson dispersion. Once the d- 
and /-bond bosons are no longer condensed, the coher- 
ence of the 6-bond bosons is built up via processes de- 
scribed by the terms bibjS ia h}j S {b\b^Si a hj S ) in the Hamil- 
tonian. 

Let us now consider the physics sketched above in 
terms of a mean-field analysis of the BFM. After perform- 
ing the corresponding decoupling procedure, the Hamil- 
tonian can be separated into three main parts: 

H = H t + H ex + H , (10) 

where H t is the hopping term, H ex the local pair ex- 
change, and Hq a term including the contributions com- 
ing from bond bosons in the chemical potential term 
and the constraint. Here, the constraint is treated 
in an approximative way, where the corresponding La- 
grange multiplier is taken as spatially homogeneous and 
which has to be determined variationally as a saddle 
point solution of the total free energy. In our mean- 
field approximation, we rewrite the first part of Eq. 
in a compact way by introducing the vector operator 
Dk = (sfc,fi sLfc 1 1 h>k,h h-k i) ano - ^s hermitian conju- 
gate D\ = (s\.p S-k,l> h\ tj h—k,l) m terms of the Fourier 
representation for the bond fermions. The single-particle 
excitations are then determined by the following mean 
field Hamiltonian: 

H t = J2 D tUkD k + c, (11) 

k 

where c = 2A— 4/i is a constant which arises in the process 
of ordering of the various bond fermionic operators. The 
4x4 matrix Uk is given by: 



"7^: 




it 


it 




-71 


it 


-it 


7 f 


it 


Ik 


A£ 


.it 


-it 


A£ 


-Ik . 



6k = — t exp[z/c • S] is the free particle spectrum of 
the original fermions, with S designating the lattice vec- 
tors linking nearest neighbor sites. 7| = €fe[0.5 * (b + 
a) 2 - f 2 } - 3^ + A, 7 £ = e k [d 2 - 0.5 * (b - a) 2 } - (j, + A 
are the single particle energy spectra for the s and h 
fermions. 7 f = rfe{[(b + a)d] - [(a - b) /]} and j sh = 

e k[(f d) + \ (a 2 — 6 2 )] are the particle/hole and parti- 
cle/particle hybridization factors between the s- and h- 
bond fermions. The pairing amplitudes for the s- and 
/i-bond fermions are A| = ek>/2[f(b + a)] and A£ = 
e k V^[d{a — b)], respectively. 

Let us next investigate how the effective hopping of 
the s, /i-bond fermions depends on the amplitude and 



the density of the condensed bosons. As one can see 
by inspection of the matrix Uk, the processes linked to 
the fermionic degrees of freedom are strongly renormal- 
ized by the strength of the condensed 6-bond bosons. 
Concerning the dispersion of the s, ft,-bond fermions, one 
can see that the contribution of the density of the d, /- 
bond bosons counteracts that of the &-bond bosons. The 
strength of the pairing amplitude being proportional to 
the product of the &-bond and d, /-bond boson conden- 
sate amplitudes, indicates the need of having both types 
of bond bosons condensed in a state where the s, ft.-bond 
fermions are paired up. The hybridization between the 
s- and /i-bond fermions manifests itself in both, the par- 
ticle/particle and in the particle/hole channel. While the 
amplitude 7^ of the s — h hybridization is linearly linked 
to the 6-bond and d, /-bond boson condensate, the pro- 
cesses related to the particle/particle mixture depend on 
the effective density of the &-bond bosons as well as the 
d, /-bond bosons. All those processes compete with each 
other, resulting in either a superconducting or an insu- 
lating state, depending on the value of g/D. 

Let us next examine the exchange and the local contri- 
butions of the BFM, which after a corresponding mean- 
field decoupling are given by: 

H ex = Lg(a 2 -b 2 ) 
H Q = -L^[4/ 2 + 2(6 2 + a 2 )] + 

L\[f 2 + b 2 + a 2 +J 2 -1} (12) 

and where L indicates the total number of bonds. The 
procedure of the present analysis is to look for a sad- 
dle point solution of the total free energy with respect 
to the bond boson amplitudes and the Lagrange multi- 
plier for the constraint. The free energy is obtained after 
performing a Bogoliubov rotation of the s, /i-bond boson 
operators in the H t term, bringing it into the diagonal 
form: 

H t ^^2E a , k vl k 7] a , k . (13) 

fc 

E a ,k (a = L--,4) are the eigen-energies of the various 
Bogoliubov quasi-particles ry a k (r) at k), obtained by diag- 

onalizing the matrix U k - They result from a unitary 
transformation of the original fermions and thus have 
the same commutations relations as Dl{D k ). The free 
energy is then given by the following expression: 

F = E M2cosh(^, fc )] +J E 

^ k i=l,..4 

E Q = g(a 2 - b 2 ) - M [4f + 2(6 2 + a 2 )] + 

X[P + b 2 + a 2 +S 2 - 1] (14) 

with (3 = l/(k B T). 

The saddle point solutions for the bond bosons to- 
gether with the Lagrange multiplier are obtained vari- 
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ationally by requiring the following extremal conditions: 
OF 

dF _dF _dF _ OF 
da db dd df 











(15) 



together with the chemical potential being fixed via the 
relation n tot = -fjj- 

In Fig|2 we report the evolution of the amplitude 
for the condensed o-bond bosons at zero temperature 
as a function of g/D. Throughout this paper we use 
a 2-dimensional fermionic tight binding spectrum efe = 
— (D/2)(cos[k x ] + cos[k y ]) for the original fermions. The 
self-consistent solutions of the non-linear Eqs. 1151 yield 
a zero amplitude for the a-bond bosons, as anticipated 
above, and an equal amplitude for the d- and /-bond 
bosons. The latter is a consequence of the particle-hole 
symmetric case considered here. We also find a clear 
separation between two distinct regions as the exchange 
coupling is varied. Below a critical g cr it/D ~ 0.4 the sys- 
tem allows a finite condensation of the d, /-bond and the 
o-bond bosons, which signifies a superconducting state. 
Yet, upon increasing the exchange coupling, the ampli- 
tude of the d, /-bond bosons gets reduced, until above 
g C riu the ground state does not contain anymore any d- 
bond bosons. The transition is continuous, merely show- 
ing a change of slope in the behavior of b. 

A particular feature of this mean field analysis at T = 
is that the order parameter d, controlling the supercon- 
ducting phase, has a finite value as soon as the boson- 
fcrmion exchange mechanism is switched on, i.e., for 
g/D ^ however small it might be. In the limiting case 
were the exchange coupling is zero, the superconduct- 
ing solution disappears because then o-bond and a-bond 
bosons are equally populated. Such a degeneracy is re- 
moved as soon as the exchange coupling is different from 
zero. The o-bond bosons are now being able to condense 
while a-bond bosons will not. Thus, in the present mean 
field treatment the system finds the solution with a fi- 
nite amplitude for the o-bond and d, /-bond bosons as 
soon as the exchange is switched on at a non-zero value. 
The size of the o-bond boson amplitude is then related 
to the constraint and via that to the total filling under 
consideration. 

At this point we would like to make a few remarks on 
the weak coupling regime. Since the mean-field analysis 
is based on a local ansatz for the saddle point solution, we 
do expect that the present approach is best suited for the 
strong-coupling limit. In the weak coupling limit such a 
scheme of approximation is not completely satisfactory. 
In that regime, the dynamics of the bond bosons and 
their delocalized behavior have to be taken into account 
in order to properly consider their phase and amplitude 
fluctuations and the related to it consequences on the 
fermionic subsystem. The bosonic fluctuations should 
renormalize down to zero the value of the critical temper- 
ature in this small coupling region. The finite amplitude 
of T c at g/D ~ 0, obtained in our present mean- field 
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FIG. 2: Phase diagram at zero temperature as given by the 
evolution of the bonding bosons and bond hole amplitude. 
SC and I stands for a superconducting and insulating state 
respectively. 



analysis, indicates that the pairing formation is driven 
by kinetics and an overestimation of the phase locking 
with respect to the amplitude fluctuations of the bond 
boson order parameter occurs. This physically unsatis- 
factory feature is corrected by including the competition 
between the phase and amplitude dynamics, as we have 
done in a recent work where an effective phase action has 
been derived after integrating out the fermionic degrees 
of freedom (see insert in Fig. 3)3i We do believe that, 
even in the weak coupling regime, a BCS type of super- 
conductivity should not occur since the charge dynamics 
intrinsically involves the two types of s— and h— bond 
fermions. A direct decoupling of the exchange interaction 
in the starting representation of the BFM would hence 
not be justified as the hard core bosons are localized fields 
and thus in principle cannot have a phase coherent am- 
plitude in the transverse mode. In an attempt to include 
the mobility of the bosons, in a self-consistent perturba- 
tive (weak-coupling limit) approach, it has been shown 
that already above the superconducting transition tem- 
perature one gets, apart from the regular free band spec- 
trum, two rather dispersion-less branches^ As we will 
see below, such an excitation spectrum resembles that 
obtained in the present approximation. Hence, within 
our saddle point solution qualitatively correct features 
can be recovered even in the weak-coupling regime. 

Concerning the dynamics of the bosons, we should 
make a few remarks about the possible nature and the 
character of the bosonic excitations near the quantum 
phase transition between the superconducting and the 
insulating state - in a frame that is beyond our present 
approximative mean field scheme of investigation. Focus- 
ing on the dynamics on the insulating side, we observe 
that the fermionic excitations are gapped, while the b- 
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FIG. 3: Phase diagram at finite temperature as a function 
of g/D as given by the evolution of the 6-bond boson and 
d-bond boson amplitude. SC, BP and N stands for a super- 
conducting, bond coherent boson phase, and normal state, 
respectively. The insert shows the evolution of the critical su- 
perconducting temperature as obtained within a path-integral 
approach p— 

bond bosons may have low energy gapless excitations 
resulting from their itinerancy. However, such excita- 
tions are not "charged" ', i.e., capable of creating a pair 
of fermions, or a localized boson, since that would require 
at the same time the presence of both 6-bond as well as 
d, /-bond bosons. Hence, two scenarios are possible as 
g/D is increased: i) the d, /-bond boson coherence and 
their density go to zero simultaneously, ii) the phase co- 
herence of the e?-bond bosons drops to zero before their 
density vanishes. Case i) corresponds to a direct transi- 
tion between the superconducting to an insulating con- 
figuration with coherent 6-bond bosons. Case ii) allows 
to have an intermediate state, where a finite density of 
d, /-bond bosons occurs in a state that is not condensed. 
In this case, the current can be carried directly by the 
mobile 6-bond and d, /-bond bosons, and the presence 
of these low energy excitations can induce a metallic-like 
state. This regime is quite exotic, since the effect of the 
a-bond bosons in activating itinerant pairs can induce a 
gap in the bosonic sector and, in this way, result again in 
an insulating state. The picture that finally could emerge 
is one of a system that undergoes first a transition from 
a superconductor to a metallic-like phase (built out of 
a state where incoherent itinerant 6-bond and d, /-bond 
bosons are present and which ultimately transits to an 
insulating state where, (driven by the exchange coupling 
and due to the constraint) there exist only condensed 
6-bond bosons and no condensed d, /-bond bosons. 

Let us now analyze the behavior of the two coherent 
states as a function of temperature (see Fig. In the 
structure of the phase diagram one notices three regions 
that are reminiscent of the configurations found at zero 



temperature. The superconducting state (SC) is charac- 
terized by a non-zero amplitude for the 6- as well as d, f- 
bond bosons, the BP (bond phase) by condensed 6-bond 
bosons but a zero amplitude of the d, /-bond bosons, 
while in the normal (N) high temperature phase none 
of the bond bosons are condensed. As a consequence of 
the competition between the 6-bond and d, /-bond boson 
condensation, one has two regimes depending on the ex- 
change coupling in the region where SC is stable at zero 
temperature. For values of the exchange coupling below 
g* / D = 0.125, the SC critical temperature (Im) and that 
associated with the 6-bond boson condensation are 
very close in amplitude so that there occurs an almost 
direct transition from the SC to the N phase. In this 
limit, the scales of energy that controls the 6-bond and 
d, /-bond dynamics are comparable. Going to a larger 
exchange couplings [g > g*), the characteristic temper- 
ature marking the onset of SC is reduced due to the in- 
creased population of the 6-bond bosons. But contrary to 
that, Tb grows almost linearly with the coupling. Thus, 
outside the quantum phase transition region there occurs 
a large region in parameter space, where the system is not 
superconducting but is still characterized by a coherent 
state of 6-bond bosons. 



V. EVOLUTION OF THE CHARGE GAP AND 
OF THE SPECTRAL FUNCTION FROM THE 
SUPERCONDUCTING TO THE INSULATING 
REGION 

In this section we analyze the excitation spectra both 
for the s- and ft,-bond fermions and for the original c 
fermions. We shall focus on the evolution of the spectral 
function when moving from the superconducting to the 
insulating region in the phase diagram. In the previous 
section, we have constructed the vector D^, with D/- } i de- 
noting its i—th component. Due to its structure, the time 
dependent correlation function for the vector D)~ contains 
the information for the matrix Green's function both for 
the diagonal and the off-diagonal anomalous part of the 
s- and /i-bond bosons. Using the relation which connects 
the s- and /i-bond-fermions and the 6-bond bosons with 
the original c type fermions, we are able to extract the 
spectral function for the original fermions. For that, it 
is convenient to introduce the following time dependent 
correlation function: 

<%(T) = -{T(p kti (r)Dl tj (pj)) (16) 

where i,j — 1, ..,4 and r is the Matsubara imaginary 
time, T the usual time ordering operator, and (...) in- 
dicates the average at finite temperature. Due to the 
bilinear structure of the Hamiltonian, we readily deter- 
mine the Green function via its equation of motion. We 
thus obtain the following expression in matrix notation: 

G fc K) = (iw n l-U k )- 1 . (17) 

The Green's function for the original c operators, is con- 
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FIG. 4: Zero temperature evolution of the charge gap asso- 
ciated with the excitation of a c fermion as a function of the 
exchange coupling. 



sequently a linear combination of the contributions of 
the different components GV (r) which are weighted by 
the various amplitudes of the condensed 6-bond bosons. 
This is because of our mean field approach where in the 
expression for the fcrmionic c operators, Eq. (J5J, the 
various boson operators are replaced by their mean field 
averages, i.e., 

- 1 - p - 

4a = PafSiB + ~^4a ( b + S ) + ~j= (« " b ) Ka + h\ a d 

Cia = Pal 4a + ^ s w (f> + 5 ) + 7^= ( 5 ~ *>) h ia + ^iad 

(18) 

Expanding the expression for the time dependent 
Green's function of the c fermions G£ CT (w„), we derive 
the following structure for it: 





G ka{u n ) = Tr[F ■ Gfc(w n )] (19) 

where Tr indicates the trace of the matrix indices, and 
F is a 4x4 matrix whose coefficients depend on the am- 
plitude of the boson condensates. 

In our investigation of the various phases of the BFM, 
let us start with by examining the evolution of the charge 
gap, as extracted from the density of states, N c (uj) — 
fdeN(e)A%(oj), with A%(w) = {-l/ir)Im[G%(w)] and 
N(e) the bare density of states. Both in the SC- and 
the I-phase the 6-bond bosons are condensed. Thus, the 
charge gap (A c ) represents an estimate of the energy re- 
quired to break a 6-bond boson and create a pair of h- 
and s-bond fermions in case those bond fermions sim- 
ply hybridize or pair up. As one can see from FigQ] 
A c has a non-monotonous behavior as the exchange cou- 
pling g is varied from weak to strong coupling. For gj D 




FIG. 5: (Color online) Dispersion of the Bogoliubov branches 
Ei and the related spectral weight at zero temperature for 
the c-type excitation as a function of the exchange coupling. 
From top to the bottom the exchange assumes the following 
values: g/D = 0.1, 0.3, 0.4, 0.6, respectively. The points of the 
two dimensional Brillouin zone are: V = [0,0], Z — [n, 0], X = 
[ir,ir]. 



9 



smaller than about 0.15, the energy to excite a single 
c-type fermion grows slightly as a consequence of the 
concomitant effect of the constraint, the hybridization 
between the s- and /i-bond fermions as well as the in- 
creasing pairing between them. Beyond g/D ~ 0.15, the 
bond-hole boson amplitude diminishes, which drives the 
charge gap to a lower bound. Upon further increasing 
g/D, the systems transits from a superconducting to an 
insulating phase. The non-monotonic behavior occurring 
in the superconducting regime is related to the compe- 
tition between the two possible channels for breaking a 
6-bond boson (see Fig. lb). Moving into the insulating 
side of the phase diagram, we notice that the charge gap 
increases almost linearly with g. This indicates that now, 
for breaking one 6-bond boson in order to create a pair 
of s, /i-bond fermions, one has to overcome the strength 
of the local exchange. In this limit, it is the constraint 
that determines the local scale of exchange energy in the 
excitation spectrum. 

Up to now we have concentrated our attention on the 
amplitude of the 6-bond bosons and their, related to it, 
charge gap and its evolution from the SC to the I phase. 
We next shall investigate how the dispersion and the re- 
lated spectral weight associated with the c-type excita- 
tions of the Bogoliubov branches get modified by tuning 
the exchange coupling. In Fig. [S] we report the behav- 
ior of the spectral function for the original fermions as 
a function of their momentum, when moving along the 
major symmetry lines of the two-dimensional Brillouin 
zone. 

We observe that the structure of the excitation spec- 
tra is qualitative different in the superconducting phase 
and the insulating one. We obtain four branches in the 
SC region and only two in the insulating one. Concern- 
ing the dispersion and the position where the charge gap 
opens, there are distinct differences in the two parts of the 
phase diagram. Starting from the weak coupling regime, 
one notices that the gap opens up at the Z point of the 
Brillouin zone. This is due to a crossing at this point of 
the fermion-like and fermion-hole-like tight-binding spec- 
trum, for the particle-hole symmetric half filled band 
case, studied here. Looking at the evolution of the spec- 
tral function, the excitations with the dominant weight 
follow a dispersion that is reminiscent of a free tight bind- 
ing spectrum, except close to the Z point, where the gap 
opens. Furthermore, in this region of coupling, two al- 
most dispersionless branches with small spectral weights 
are present. They arise from the hybridization between 
the s- and /i-bond fermions. Their small spectral weight 
is a consequence of a counter-active action between the 
6-bond and d, /-bond bosons in renormalizing the hy- 
bridization between the bond fermions. Increasing the 
exchange coupling induces a modification in the spectral 
weight distribution and in the width of the dispersion. 
The dominant Bogoliubov branches get more dispersive 
upon approaching the quantum phase transition due to 
an increase in the 6-bond boson amplitude 6. Their spec- 
tral weight however diminishes down to zero due to the 



reduction of the bond-hole and double paired boson am- 
plitudes d and / as we approach the quantum critical 
point. This change is accompanied by a partial redis- 
tribution of the spectral weight in the lower and upper 
Bogoliubov branches. Crossing the critical exchange cou- 
pling, the bands below (above) the chemical potential 
develop a unique structure and one no longer discerns 
any trace of the four Bogoliubov branches which were 
caused by the pairing and the hybridization between the 
bond- fermions. In this region, only the process of mix- 
ing the s- and /i-bond fermions is active. Now, a 6-bond 
boson can only break into a s- and ft.-bond fermion but 
the energy cost for that grows proportionally to the ex- 
change coupling and is controlled by the constraint. The 
emerging bond-fermions have propagating features and 
hybridize with each other which results in a dispersion 
that is renormalized by the 6-bond boson amplitude 6. 



VI. CONCLUSIONS 

In the present study we have illustrated the evolution 
of pairing correlations in a two-component system where 
local bond pairing is induced in an ensemble of itiner- 
ant fermions. We have shown that there are two possible 
bond pair configurations which can materialize. Depend- 
ing on whether we have a finite or a zero amplitude of 
the bond- hole (respectively double-paired-bond) bosons, 
we obtain a superconducting or an insulating phase. The 
bond fermions which emerge in such system control the 
occurrence of the one or the other bond pair liquid. Those 
bond fermions are both: (i) single fermions on a bond and 
(ii) fermions attached locally to a boson on that same 
bond. The dynamics of these bond fermions and the in- 
terplay between the related pairing and the hybridization 
process crucially determines the competition between the 
superconducting and the insulating paired liquid states. 

To extract the main features of this scenario, we have 
used a bond operator formalism, which has been widely 
applied in the literature for the spin analogue quantum 
problems. In the present problem, it has turned out to be 
particularly useful because of the intrinsic nature of the 
dimer formation of resonant pairs. Moreover, though the 
bond operator method, within the present mean field ap- 
proach, is best suited for the strong coupling limit, it can 
indeed describe the quantum phase transition because, as 
we have shown, it occurs in a regime of coupling where 
g/t ~ 2. This, a posteriori, justifies the use of an ap- 
proximation that is based on a local saddle point ansatz 
for the bond bosons. Still, as we have discussed in the 
Introduction, it is already within the local structure of 
the bond configurations that one recovers the objects to 
distinguish between the superonducting and the insulat- 
ing state. Among the non-bonding configurations, there 
is the bond-hole state, which turns out to be the relevant 
bond boson, because the formation of the superconduct- 
ing state is attributed to the simultaneous presence of a 
phase coherent amplitude for the bonding bosons as well 
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as for the bond holes. This is a key element that allows 
to follow the evolution from the insulating (condensation 
of the pure bonding bosons) to the superconducting state 
(pairing of the s and h fermions, due to the phase coher- 
ence of the bonding bosons and the bond holes). 

We have obtained the phase diagram as a function of 
the boson-fermion exchange coupling and by varying the 
temperature. We have shown how, with increasing the 
exchange coupling, the continuous reduction of the am- 
plitude of the bond-hole condensate drives the quantum 
phase transition between the superconducting and the 
insulating state. This behavior is dictated by the inter- 
play between the dynamics of the bond fermions and the 
constraint controlling the occupation of the various bond 
operators. In this first attempt to capture the salient 
aspects of that scenario, on the basis of a mean field ap- 
proach in the bond operator formulation of it, we have fo- 
cused our attention on the excitation spectra and the na- 
ture of the charge gap. In the superconducting region the 
single particle gap exhibits a non-monotonic evolution 
as the exchange coupling is increased and the quantum 
phase transition is approached. This is a sign of the dou- 
ble nature of the dynamical variables at play and which 
contains features of pairing as well as of hybridization 
of those bond fermions. A completely different behavior 
characterizes the insulating regime where it is the local 
exchange which sets the scale of energy. Coming from 
the insulating side, entering the superconducting phase 
one is led to the picture where the insulating charge gap 
continues to exist in the superconducting phase as one 
of the components of the charge gap, and whose ampli- 



tude diminishes as g goes to zero. This could signify a 
breakdown in cascade of the order in such systems, as g 
is increased. First destroying the superconducting phase 
but without the amplitude of the pairs up to the quan- 
tum critical point and then destroying gradually the pair 
amplitude when g is sufficiently large in the insulating 
state. 

Concerning the excitation spectra, there are four Bo- 
goliubov branches in the superconducting region and 
which continuously evolve into a two branch configura- 
tion in the insulating region. For a 2-dimensional system 
and the corresponding tight-binding spectrum for the 
original bare fermions, the pairing gap opens up at the 
Z point of the Brillouin zone. The interplay between the 
pairing and the hybridization between the bond fermions 
also controls the spectral weight redistribution between 
the various Bogoliubov branches as the exchange cou- 
pling is tuned through the quantum phase transition. 

Further analysis on the issues which have been dis- 
cussed here is presently in progress and deals with such 
questions as the effects expected for hole doping away 
from the half-filled band case and of the feedback on the 
bond fermions arising from the dynamics of the bond 
bosons. A major attempt in this direction will concern 
the case when the system is in the proximity of the break- 
down of the coherent superconducting state where possi- 
bly a metallic phase of bosonic fermion pairs could exist. 
For that purpose the present bond-mean-field analysis 
will be improved by taking into account the phase fluctu- 
ations of the bosonic mean field coherent states desribed 
here. 
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